Bounding the set of quantum correlations 
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We introduce a hierarchy of conditions necessarily satisfied by any distribution P a p representing 
the probabilities for two separate observers to obtain outcomes a and (3 when making local measure- 
ments on a shared quantum state. Each condition in this hierarchy is formulated as a semidefinite 
program. Among other applications, our approach can be used to obtain upper-bounds on the quan- 
tum violation of an arbitrary Bell inequality. It yields, for instance, tight bounds for the violations 
of the Collins et al. inequalities. 



The correlations between two separated physical sys- 
tems can be characterized by the joint probabilities P a p 
that an observer who performs a measurement X on the 
first system gets an outcome a £ X and that an observer 
making a measurement Y on the second system gets an 
outcome j3 € Y. If the observed system is in an en- 
tangled state, these joint probabilities may violate a Bell 
inequality, implying that quantum theory is not, in Bell's 
terminology, locally causal [H. Although quantum cor- 
relations are not constrained by Bell's locality principle, 
they are not arbitrary since a general joint distribution 
P a p cannot always be viewed as originating from mea- 
surements performed on a shared quantum system Q ■ 

In this paper, we investigate the restrictions on bipar- 
tite correlations imposed by the quantum formalism. The 
question that we seek to answer is the following: Given 
an arbitrary distribution P a p, do there exist a quantum 
state p on a joint Hilbert space Ha ®'Hb and local mea- 
surement operators E a = E a ® I and Ep = I ® Ep, such 
that P aP = tr {E a Ep p)l 

From a fundamental point of view, one motivation for 
studying this problem is simply to understand which cor- 
relations can arise between two systems within our cur- 
rent description of nature. Another is to develop tools 
to detect the possible non-quantumness of some set of 
empirically observed correlations. Answering the above 
question is also of practical interest for various applica- 
tions in quantum science, for instance, for the design of 
nonlocality tests more resistant to experimental noise and 
detector inefficiencies. In general, characterizing the set 
of quantum correlations is essential to understand better 
the extent to which quantum mechanics is useful in infor- 
mation processing tasks such as communication complex- 
ity Q and key distribution 0, [f| . An usual problem in 
these contexts is to determine what is the maximal viola- 
tion of a Bell inequality allowed by quantum mechanics. 

Answering the above question is in general a difficult 
task; the simple strategy of searching over all quantum 
states p and measurement operators E^, which in prin- 
ciple can be of arbitrary dimension, is clearly unfeasible. 
The first to tackle this problem was Tsirelson who 
derived the maximal quantum violation for the Clauser- 
Horne-Shimony-Holt (CHSH) inequality Q. Following 



his seminal work, many attempts have been made to un- 
derstand better how the quantum formalism constrains 
the correlations between two parties Q , but most of the 
obtained results only apply to a restricted set of situa- 
tions or lack an efficient generalization. Among these, we 
mention the works of Landau 0] and Wehner [lj| who 
realized that the problem of deciding if a set of "correla- 
tion functions" obtained from measurements with binary 
outcomes admits a quantum representation can be cast 
as a semidefinite program (SDP). Their approach is in- 
teresting in that semidefinite programs are convex opti- 
mization problems for which powerful com put ational and 
theoretical methods have been developed 11 1 . 

We show here how to design criteria that distinguish 
correlations that can be reproduced through local mea- 
surements on a quantum state from those which can- 
not. Testing any one of these criteria amounts to solve 
a SDP. Contrary to previous constructions, our method 
can be applied to any distribution P a p of joint probabil- 
ities, corresponding to a configuration with an arbitrary 
number of measurement choices and outcomes. Seen as 
a whole, the tests that we introduce exhibit a hierarchi- 
cal structure, in the sense that they can be organized as 
a sequence of conditions, each condition in the sequence 
being stronger than the preceding one. As an illustra- 
tion of the effectiveness of our approach, we present two 
applications of it. First, we derive a non- linear inequal- 
ity satisfied by quantum mechanics which strengthens a 
previous inequality due to Tsirelson 12 1, Landau 0] and 
Masanes [T3|. As a second application, we give a tight 
bound for the violation of the Collins et al. inequali- 
ties \u\. 

Preliminaries. Before entering in the details of our 
construction, let us first give some definitions and state 
the assumptions made through this paper. We assume 
that the two parties, Alice and Bob, choose their mea- 
surements from a finite set of possibilities, and that each 
measurement may yield one out of a finite set of out- 
comes. Note that we think of outcomes corresponding 
to different measurements as being labeled distinctly, so 
that each outcome a of Alice (or j3 of Bob) is unambigu- 
ously associated to a unique measurement X (or Y). 

Refining the statement made earlier, we say that a dis- 
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tribution P a p admits a quantum representation if there 
exist a joint quantum state p on "Ha ® H b , a set of pro- 
jection operators E a = E a ® I acting on Alice's system 
and a set of projection operators Ep = I ® Ep acting on 
Bob's system, such that 



P a/3 = tr (E a E p) 



(1) 



Projectors corresponding to outcomes belonging to the 
same measurement M should (i) be orthogonal: E^Ey = 
for p, v £ M, p v, and (ii) sum to the identity: 
J2fi£M En = I- By definition, we also have that (Hi) 
E^ = Ej l = E^ and that (iv) the projectors on Alice and 
Bob's side commute with each other: \E a ,Ep] = 0. 

Note that the most general description of a quantum 
measurement corresponds to a positive operator valued 
measure (POVM) rather than a set of projection opera- 
tors. But since a POVM can be viewed as a projective 
measurement on a system of larger dimension, and since 
we do not impose any constraints on the dimension of 
the Hilbert space Ha ® Wb, no generality is lost with 
our definition. 

Necessary conditions for quantum probabilities. We 
now introduce a family of conditions satisfied by any 
quantum distribution P a p. We thus start by assuming 
that there is a quantum state p and a set {E^} of projec- 
tion operators satisfying Eq. |T]) and properties (i)-(iv), 
and seek new implications from these assumptions. 

By taking products of the operators E^ and linear 
combinations of such products, we can define new op- 
erators, for instance E a EpE a > or J2 a c aE a (note that 
these new operators are not necessarily projection op- 
erators anymore, nor even hermitian operators). Let 
S = {Si, . . . , S n } be a set of n such operators. Asso- 
ciate to the set S a n x n matrix T through 



Tij = tr [SjSjp 



(2) 



By construction, the matrix T is hermitian, it satisfies 
the linear identities 



h3 



if ^ CijSjSj = X d a pE a Ep , 

a.,/3 



(3) 



(4) 



and it is positive semidefinite, 

r y o. 



(5) 



The linear constraints ([3]) directly follow from the linear- 
ity of the trace in ([2]). The important point is that they 
partly reflect the properties (i)-(iv) satisfied by the op- 
erators E^. For instance suppose that S contains an op- 
erator Sj = E^ and a subset of operators {Sk | k £ K,} = 



{E^Ey | v £ M} for some measurement M. Then, prop- 
erty (ii) implies T,kefC S k = T, v eM E n E » = E » = S j 
and thus J^keic = • ^ s another example, sup- 
pose that Slj = E a and Sj = EpE a i with a,a' £ X and 
a 7^ a'. Then, using successively properties (Hi), (iv) 
and (i), we find SjSj = E a EpE a > = E a E a iEp = 0, and 
thus = 0. The conditions ((4|) arc obtained by making 
use of ([1]) in ([2|) and relate the entries of the matrix T to 
the specific set of probabilities P a p under consideration. 
Finally, to establish (J3J), remember that an n x n matrix 
r is semidefinite positive if and only if v'Tv > for all 
v £ C". Expanding this expression, we get 



v f Tv = v*tr [SjSjpj Vj 



>o, 



(6) 



since p is a positive operator. 

We thus just showed that for any quantum distribu- 
tion P a p , there necessarily exists for each set S a matrix 
r satisfying the linear constraints Q and ^ and the 
condition of positivity If for some iS it is not possi- 
ble to find a matrix T satisfying these properties, we can 
therefore conclude that the correlations characterized by 
the distribution P a p cannot be reproduced through local 
measurements on a quantum state. Determining if there 
exists a matrix T satisfying conditions (J3j> , j4|, and (J5]) 
amounts to find a positive semidefinite matrix satisfying 
a set of linear constraints and is a typical instance of 
semidefinite programming (see 11( for a review of SDP). 
All the techniques developed in this context can thus be 
applied to evaluate our conditions. 

To give a concrete example of our method, consider 
the case where S = {E a } U {-Eg} is simply the set of 
projectors of Alice and Bob. Suppose that there are m 
different measurement outcomes a = 1 , . . . , m for Alice 
and m different outcomes = m + 1, . . . , 2m for Bob. 
Then S = {E\, . . . , E m ,E m +i, ■ E 2m } and applying 
the construction defined by ([2]) , we find that T is a 2m x 
2m matrix of the form 



Q P 
P T R 



(7) 



where the submatrix P is simply the m x m table of 
probabilities with entries P a p, and the submatrices Q 
and R satisfy 



}aa — Pa, 



R 



Pn 



Qaa' = (a, a £ X and a ^ a') , 
Rpp, = (13,13' £Y and 13^(3'), 



where P a = J2peY P a/3 and Pp = J2 a ex P <*i3 are the 
marginal probabilities for Alice and Bob, respectively. 
The form of the matrix {7J is defined by the linear con- 
straints Q and ([?]). The only entries of T which are 
not determined by these constraints are the entries Q aa ' 
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with a G X and a' G X' belonging to different mea- 
surements of Alice (X 7^ X'), and the entries Rpp with 
/3 € y and /?' 6 F' belonging to different measurements 
of Bob (F ^ F' ). In a quantum scenario, these entries 
would correspond to non-commuting measurements per- 
formed on each subsystem and would thus be unobserv- 
able. Nonetheless, if the correlations P a p have a quan- 
tum origin it is possible to assign values to these undeter- 
mined entries such that the overall matrix |7]) is positive 
semidcfinite, in accordance with ([5|). As mentioned ear- 
lier, semidcfinite programming can be used to determine 
if the matrix (7} can be completed in such a way. 

A hierarchy of conditions. We have shown how to de- 
sign tests that distinguish quantum from non-quantum 
correlations. Each set S of operators used in our con- 
struction yields a different condition, and the choice of 
a particular S may thus seem arbitrary. But in fact not 
all conditions built according to our instructions are in- 
dependent. Moreover, they can be organized in a hierar- 
chical structure, such that they can all be checked in a 
systematic way. 

To see this point, first note the easily established fact 
that if every operator in a set S can be written as a linear 
combination of operators in another set iS', then the con- 
ditions obtained from S' are at least as constraining as 
the one obtained from S. By this we mean that if there 
exists a matrix T' satisfying the constraints ((3j) , ((4|), (J5]) 
associated to S', then there also exists a matrix V satis- 
fying the corresponding constraints associated to S. The 
second observation is that the set T m = {E^ 1 . . . E^ m } of 
all possible products of m projectors generates by linear 
combinations all the operators that are linear combina- 
tions of products of ml projectors, with ml < m. 

Using these two properties, it is possible to remove 
the apparent arbitrariness on the choice of the set S, 
and check all the conditions that can be built using our 
method in a comprehensive way. Start first with the con- 
dition based on the set T\ = {E a }U {Ep}, which consists 
only of the projectors of Alice and Bob, and which we 
denote shortly as %_ = If this condition is satis- 

fied, consider the bigger set T 2 = {E^E U } consisting of 
all products of two projectors. If the corresponding test 
is passed, move to T 3 = {E^E^Er}, and so on. Since the 
elements of T n -\ can be generated by linear combinations 
of the elements of T n , each successive test is at least as 
good as the previous ones. The sets T n , then, define a 
hierarchy of necessary conditions for testing the quantum 
origin of P a p- Note that if a test fails at some point in 
this hierarchy, we can immediately conclude that the cor- 
relations under consideration are not-quantum and there 
is no need to proceed with the successive tests. The con- 
dition based on the matrix ([7]) corresponds to the first 
test in this infinite hierarchy As an illustration, we now 
present two applications of our method. 

Application 1. The first example involves two mea- 
surements for Alice, X = 1, 2, and two for Bob, Y = 3,4, 



where each measurement yields one out of two outcomes, 
+1 or — 1. This situation is thus characterized by six- 
teen probabilities P(±x)(±Y), to which we can associate 
eight projectors E±m (M — 1, ... ,4). Alternatively, one 
can characterize this scenario by specifying the correla- 
tion functions Cxy = J2 a b a b P{aX)(bY) and the marginal 
quantities Cx = J2 a a ^ aX aim = J2b ^ ^>Y for each 
measurement. The first test in the hierarchy corresponds 
to the set T\ = {E±\, . . . , E± 4 }, or, equivalently, to the 
set S = {I, a\, . . . , 04}, where a M = E+m — E-m, which 
is linearly equivalent to T\ . The corresponding 5x5 ma- 
trix r is 



(8) 



where we have only given the upper triangular part 
of r since it is hcrmitian. The parameters u, v corre- 
spond to entries that are not determined by our construc- 
tion; but if the correlations represented by the quantities 
{Cx,Cy,Cxy} are quantum, it is possible to find val- 
ues for u and v such that the matrix ((8|) is semidcfinite 
positive. 

Note that a necessary condition for the matrix §E§ to 
be semidefinite positive is that the bottom-right 4x4 
submatrix satisfies 



/I Ci 


c 2 


c 3 


cA 


1 


u 


C13 


C14 




1 


C23 


C*24 






1 


V 


V 






1 / 



/I 



V 



U C13 Ci4^ 
1 C23 C24 
1 V 



>- 0. 



(9) 



1 / 



This condition for the correlation functions C xy w as al- 
ready introduced by Landau and Wehner |lol | . It is 
shown in Q that there are values u and v such that © 
holds if and only if the correlations Cxy satisfy the in- 
equality 

|asinCi3 + asinCi4 + asinC*23 — asinC24| < it , (10) 

and the three other ones obtained by permutation of 
the measurements. These inequalities have also been de- 
rived from a different perspective by Tsirelson [l2T ] and 



Masanes [13 1. 



The correlation functions Cxy, however, only provide 
partial information about the full probability distribu- 
tion. Our construction ([8|) , on the other hand, deals with 
the full probability distribution, including the marginals 
Cx and Cy . Using the same techniques as used in [9( to 
derive (I10p from ([9]), it can be shown that the condition 
that the matrix JHJ) is semidcfinite positive is equivalent 
to the inequality 

|asinDi3 + asinZ?i4 + asinl?23 — asinL^I < 7r (11) 



4 



where D, 



(di - CiCMyfil-CfKl-C]), 



and to 



the inequalities obtained from by permutation of 
the measurement choices. If we neglect the marginals 
by imposing Cm — we recover (fl"0|) . As a test on the 
full distribution, however, our inequality is stronger than 
(fT0]l since there are probability distributions that satisfy 
(fTH)) but which violate (fTTj) . Yet, (fTT|) is not a sufficient 
condition for a full probability distribution to admit a 
quantum representation, as we have examples of correla- 
tions that satisfy (fTTj) but which fail the successive step 
in the hierarchy. 

Application 2. By maximizing the violation of a Bell 
inequality over the set of probability distributions satis- 
fying one of our conditions, we obtain an upper-bound 
on the violation of this inequality by quantum mechan- 
ics (since such conditions are satisfied by every quantum 
distributions). A Bell inequality is a linear combination 
of the probabilities P a /3, and since these probabilities are 
related in a linear way to the entries of the matrices F, 
obtaining such an upper-bound can be cast as a SDP. 
Let us stress that this bound is an upper-bound on the 
global maximum of the inequality since we do not sup- 
pose anything about the measurements and states except 
that they are quantum. 

As an example, consider the CHSH inequality, which 
reads C13 + C14 + C23 — C24 < 2. Maximizing this ex- 
pression over all distribution satisfying the criterion ([9]) 
presented in Application 1 corresponds to the SDP, 



maximize C13 + C14 

subject to r y . 



C23 — C-2 



(12) 



The solution to this optimization problem is 2\/2, as 
noted by Wchner [lo| . and we thus recover the well- 
known Tsirelson bound. More generally, for any given 
Bell inequality, SDPs can be associated to each of the 
conditions T\ , T~2 , . . . of our hierarchy, the solutions of 
which would yield a sequence I\ > 1% > ... of upper- 
bounds on the quantum violation of the inequality. Note 
that a tight bound may already be reached after a finite 
number of such iterations, as the CHSH example shows. 

We have applied the approach just outlined to the 
Collins et al. inequalities [14[. This family of inequal- 
ities involves two measurement choices per party and d 
outputs per measurement, and can be viewed as a gener- 
alization of the CHSH inequality for systems of dimension 
greater than two. In [lit , lower-bounds for the maximal 
violation of the Collins et al. inequalities are given for 
d = 3, . . . , 8 by exhibiting a particular set of measure- 
ments and states of dimension d x d yielding high viola- 
tions of the inequalities. These quantum states have the 
particularity to be non-maximally entangled. For d = 3, 
the reported violation is I* = 2.9149, the local bound 
of the inequality being Ii oc < 2. We have numerically 
solved the SDP corresponding to the first tests in our 
hierarchy using the SeDuMi Matlab toolbox (l6| . For 



d = 3, the condition 7i yields the bound I < 3.1547, 
which is about 10% higher than the violation reported in 
[Tit - The second condition T 2 , however, yields the bound 
I < 7„ = 2.9149, proving that the partially entangled 
state and the measurements described in J15J are opti- 
mal. We have also solved the SDP for d = 4, . . . , 8. As 
for the d = 3 case, the first tests in the hierarchy are 
about 10% above the values of [l5JI, but the second tests 



give the same results as the ones reported in [15J , demon- 
strating that these are the optimal quantum violations. 

Conclusion. The approach given in this work opens 
a new way to study the correlations between two sepa- 
rate quantum systems. There are several possible exten- 
sions of our technique, for instance to systems of more 
than two parties, and many potential applications of it, 
among others to study non-local properties of quantum 
correlations, such as their monogamous character [l7T |. 
Two questions that remain open is whether the hierar- 
chy of conditions that we have introduced is complete, in 
the sense that a set of correlations satisfies every condi- 
tion in the hierarchy if and only if it admits a quantum 
representation, and, if this is the case, whether it is in 
general necessary to check all the conditions in this in- 
finite hierarchy or whether it is sufficient to stop after a 
finite number of steps. 
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